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Abstract 



Non-perturbative results are obtained for multi-point correlation functions 

of the model of (2 + l)-dimensional relativistic fermions in a random non- 

Abelian gauge potential. The results indicate that the replica symmetry for 

this model is unbroken. We calculate the diffuson propagator and show that 

DC-conductivity for this model is finite. 
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Typeset using REVTpjX 



In this letter we continue to study the model of (2 + l)-dimensional relativistic fermions 
interacting with random non-Abelian gauge potential. This model was introduced in Ref. 1 
in relation with effects of disorder in two dimensional d-wave superconductors. It was shown 
that after averaging over the static disorder one gets the following Euclidean action: 

S = S + J d 2 x[ieifj ayP ifj ayP + %U ^ ° ^ a , P A^a, 9 ] (1) 
S = f d 2 x[^ pllx d^ a , p + cJ;j;] (2) 

Spinor fields ip ajP with p — 1, ...R and p — R + 1, ...2R describe Fourier components of the 
original fermions with frequencies e + uj/2 and e — u/2 respectively, where R — > is the 
number of replicas. Greek indices run from 1 to N. The flavour currents are given by 

2R 

p=l 

with r a 's being generators of the SU(N) group. In the context of disordered superconductor 
N denotes the number of nodes of the order parameter on the Fermi surface. The coupling 
constant c is related to the random potential. The matrix A 3 is the diagonal matrix with 
R elements +1 and R elements —1 on the main diagonal: A 3 = diag{l, ...1; —1, ... — 1). In 
what follows we shall also use the matrix A 1 whose only nonzero elements are situated on 
the main antidiagonal: A 1 = antidiag(l, 1, ...1). 

The Abelian version of this model has been studied by Ludwig et. al 2 . The density of 
states (DOS) is expressed as 

p(e) = -Me lim RT l Tr < $tx)ib(x) > L=o (3) 

The main difference between the present problem and the conventional localization is that 
in the present case DOS is strongly affected by the disorder. Being linear in energy in the 
absence of disorder (p(e) ~ |e|), DOS changes its behavior at low energies p(e) ~ \e\ v where 
v = 1/(2N 2 — l) 1 . It will be shown later that with DOS vanishing at zero energy the 
low energy states remain delocalized. Another difference is that the conductivity of (2 + 
l)-dimensional massless Dirac fermions with an arbitrary dumping is constant 2,3 : 
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(4) 



The incensitivity of the conductivity to dumping makes one suspect that it may remain 
finite in the presence of disorder. We shall prove that this is indeed the case. 

The model (1,2), as it stands, contains fast and slow degrees of freedom. The latter are 
SU (iV)-singlets described by the well known Q- matrices (2R x 2i?-matrices) . The symmetry 
of the model dictates that Q-matrices belong to the U(2R) — U(l) x SU(2R) group. The 
integration over the fast degrees of freedom performed in Ref. 1 gives the effective action 
for the matrices Q = gexp^Q} (7 = J2tt/NR) in the form of the Wess-Zumino model on 



W[SU(2R); 9] = j^J d 2 xTr(d^g + d^g) + — J d 3 xe abc Tr(g + d a gg + d b gg + d c g) (5) 



The quantity M is the energy scale introduced by the disorder M ~ exp[— 2ir/Nc], which 
marks the crossover from the bare DOS p(e) ~ |e| to the renormalized DOS p(e) ~ \t\ v . In 
this Wess-Zumino model M serves as the ultraviolet cutoff. It is worth remarking that the 
Wess-Zumino action is not equivalent to the action with the topological term which describes 
the conventional localization problem in a strong magnetic field. 

The Wess-Zumino model is well studied at finite R, where it is a critical theory with 
power law decay of correlation functions. Our primary objective is to show that this model 
remains well defined in the replica limit, that is its correlation functions have finite values at 
R — > 0. The closest implication of this fact is that the replica symmetry remains unbroken. 

In determining the replica limit we follow the general principle: calculating any N- 
point correlation function F 2 r(1, 2, ...N) one treats R as an arbitrary integer number on all 
intermediate steps of calculations until the finite expression is obtained. The replica limit is 
then defined as follows: 



the group U(l) x SU(2R): 



S = S + MeTr(Q + Q + ) + iM—Tr[A(Q + Q+)] 




F(l,2,...N)=hm^F 2R (l,2,...N) 



(6) 
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Let us study the correlation functions of the Q-fields. The problem of indices is simplified 
by the fact that Q pq matrices are slow parts of the operators 



N 



Qpq ~ ^ta.p'^L.a.q (7) 



a=l 



This fact provides us with a simple recipe: any iV— point correlation function has the same 
index structure as the iV-point function of the J2 a '0H,a,pV'i,a,g-fi e lds in the theory of massless 
free fermions. The simpliest example is the 2-point function: 

< Q Piri {z, z)Q^ 2P2 {0, 0) >= 5 rir2 5 PlP2 ^ M ^ 4Al ( 8 ) 
where Ai is the conformal dimension of the composite operator Q given by the sum of the 



dimensions of the bosonic exponent exp[iy 2-k/NR§] and the operator field g pr from the 
fundamental representation of the SU(2R) group. Using the results of Ref. 4 we get: 

Al = MiM + {2R ~ l ' m)l[N + m)] = 2^ (9) 

In the replica limit we get from ( ||) 

D{z, z) = lim R- 1 < Tr[h x Q(z, z)]Tr[A 1 g + (0, 0)] >= (M\z\y 2/N2 (10) 

R — ^0 

We identify this correlation function with the propagator of diffuson in the zero frequency 
limit. Its Fourier transformation gives the propagator in momentum space: 

D(w = 0lS )~«^ (id 

where D is the bare diffusion constant. 

We can identify other operators and their conformal dimensions. Thus, the operator 
expansion of two Q matrices contain the symmetric and antisymmetric fields 0$ and Oa'. 



Q(1)Q(2) = expizJ2n/RN<5>(l)]g(l)ex P [iJ27r/RN$(2)]g(2) 



exp[z v /87r/ J RiV$(l)]O s (l) + exp[i^/87r/ 'RN$(1)]0 A (1) + ••• (12) 

where dots stand for nonsingular terms. Taking the corresponding conformal di- 
mensions from we get the conformal dimensions of the composite fields Os,a — 



exp[iJ8Tr/RN$\Os,A i n the replica limit: 



v r 1 (R-1)(2R + 1), 2-N 
Aa = Mm + R(N + 2R) ] ~ -JP~ 
.. r 1 (R+1)(2R-1). 2 + N 
= Mm + WT2R) ] = — (13) 

We see that that the scaling dimension of the antisymmetric field is negative for N > 2, 
which indicates that the replica limit of the model (1,2) is a nonunitary theory. The next 
field is the adjoint field O a d appearing in the operator expansion of QQ + : 

Q(l)Q + (2) ~ g{l)g + {2) ~ Tr[r a g + r b g} = Of d (14) 

According to Ref. 6 the corresponding conformal dimension is A a( i = 2R/(2R + N). It 
vanishes in the replica limit and, as we shall see later, this field just does not exist in the 
replica limit. 

The central charge of our theory is the sum of central charges of the free bosonic field 
(C = 1) and the Wess-Zumino-Witten model on the SU (2R) group: 

j V (4g-l) = 2j 
N + 2R N K J 

Thus the resulting central charge vanishes; however, according to the definition of the replica 
limit ( H|) the physical correlation function of the stress-energy tensors does not vanish: 

< T(z)T(0) >= lim —^-r = -^—r (16) 
v ; v ; R^o2R2z 4 2Nz 4 V ; 

Let us now study the four point correlation function of the Q fields and show that it has 
the correct replica limit. For this purpose we use the result obtained in Ref. 4: 



< Q Pl n( z l, Zi)Q^ 2p2 (z 2 , Z 2 )Qp 3 r 3 (z3, z 3)Qt 4 p 4 ( Z 4:> Z i) >~ 

2/N 2 

]ZuZ23] , (F + F), (17) 

F [^pip2^P3P4 ^91 92 ^9394^*11 ^) ^PlP4 < ^P3P2^9l 94^9392^22(^5 ^)]) (^^) 

F = [5p 1 p 2 Sp 3 p 4 S qi q 3 5q 2 q 4 Gi2(x i x) ~\~ Sp 1 p 4 5p 3 p 2 5 qi q 2 Sq 3 q 4 G2l(x i x)] (19) 



where 
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^12^34 _ ^12^34 

X = , X = - — — 

^14^32 Z\^Zj, 2 

and the functions Gab,A,B = 1,2 are known functions. Now note that in our theory we 
shall deal only with correlation functions of TrQ, TrA 3 Q and Trh}Q. Since all correlation 
functions must be proportional to R, only the term with F (that is the term with all indices 
equal) survives in the replica limit, the term with F being proportional R 2 . Therefore we 
have 

l/N 2 



R- 1 < TrA 1 Q{l)TrQ + {2)TrQ{3)TrA 1 Q + {4) >- 



I ^14^23 I 



X 



1^12^14^13^24 1 

[xF(l/N, -l/N, 1, x)F(l + l/N, 1 - l/N, 2, x) + (x -> x)} (20) 

In order to get the operator algebra of the model we consider various limits of this 
formula. In the limit z 12 = Z34 = 5 — > we get 

D(l, 1 + 0; 3, 3 + 0) = (S)- 4 / N2 + \- + i- (21) 

^13 ^13 

From this we derive the following operator expansion: 

Q(l)Q + (2) ~ \z 12 \- 2 ' N \l + z 12 J{2) + z 12 J(2) + ...) (22) 

where J, J are the left and right current operators. Notice that the primary operator in the 
adjoint representation present for R 7^ does not appear in this expansion. 
In the limit Z13 = z 2 a = 6 — > we have 

D(l, 2; 1 + 0, 2 + 0) = \zf N - 2 ^C x + \z\-^ N+2 ^ N2 C 2 (23) 

(Ci,C 2 are constants) which shows the presence of the symmetric and antisymmetric fields 
with conformal dimensions ( O) in the operator expansion of QQ and Q + Q + . These oper- 
ators describe mesoscopic fluctuations of the local density of states. 

Now let us discuss a crossover to small momenta. As follows from Eq.( ) the operators 
Tr(l ± A 3 )Q are relevant perturbations with scaling dimensions 2 — 2Ai = 2 — l/N 2 . 
Therefore at finite u one can conjecture the following form for the diffuson propagator: 
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D <=o( w > q) = -fkj&f (^=W) ( 24 ) 

Eq.( shows that the spectrum of the diffusive modes iuj = D(u)q 2 ~ g 2-1 /^ 2 and the 
diffusion coefficient scales as D That is the quantity p(uj)D(uj) 

stays constant under renormalization (the presence of logarithms or finite corrections is not 
excluded). Since from the other side the conductivity is proportional to a xx ~ p(0)D(0), 
the absence of renormalization of p{u)D(ui) at finite frequences makes it likely that the 
conductivity remains finite in the presence of weak disorder. 

We can go futher than this qualitative analysis. Namely, using the peculiar properties 
of the action ( |5|) we can calculate the propagator of diffuson in the limit q « \uj\ l l l ~ Al . 
Let us introduce the following change into the action ( [5]): 

Tr(f ± A 3 )(Q + Q+) -> Tr(f ± A 3 )(h + Q + Q + h) (25) 

where h(x) is an external source - a matrix from the U(2R) group. It turns out that in the 
limit q « |c<j| 1 / 1_Ai we can integrate over Q-variables and get the effective action for h(x) 
in the form of expansion in |o;| -1 / 1-Al . Indeed, let us make a shift: Q = h + Q,Q = hQ. 
Using the well known property of the Wess-Zumino action, we write 

S = So(hQ) + MeTr(Q + Q + ) + tM^Tr[A 3 (Q + Q+)\ = 

S(Q; e, u) + S (h) + — / d 2 xTr[h + d^hQd^Q+] (26) 

o7T J 

Since fluctuations of the Q-field are governed by the Wess-Zumino action perturbed by 
relevant operators, they have a finite correlation length £ ~ (|o;|) -1 ^ 1_Al ). Therefore the 
last term in the action ( |2B[) gives corrections to Sq(K) which contain powers of £V/i. At 
small momenta such perturbations can be neglected. Thus at small wave vectors Sq{K) 
effectively becomes a generating functional of correlation functions of Q. In particular, we 
have 

M 2 N a 2 

K(e = ^u-q) = — < Tr(l + A 3 )A 1 QTr(f - A 3 )A 1 Q + >= --^(1 + 0{q 2 /\s\ 1 -^)) 

(27) 



This means that the crossover function f(x) in Eq.( |24[) behaves as 



/(*) 



N 1 



2ir x 



at large x. This result allows us to calculate the DC-conductivity a xx . Using the definition 
given in Ref. 2 we get from ( ^) 



Comparing this result with Eq.(4) we see that the DC-conductivity remains unrenormalized 
and therefore the disorder manifest itself only at large momenta. In the limit q — > there 
are no singular contributions. This is not typical for conventional disordered systems where 
disorder changes a statistics of low-lying energy levels. In the same time one probably should 
not be surprized by our results. After all they show that the model with disorder remains 
integrable and we know that integrability and fractality of the spectrum are inconsistent. 
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